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Abstract. 

Horizontal branch (HB) stars play a particularly important role in the "age debate," since 
they are at the very center of the long-standing "second parameter" problem. In this review, I 
discuss some recent progress in our understanding of the nature and origin of HB stars. 

Keywords, stars: abundances, stars: evolution, (stars:) Hertzsprung-Russell diagram, stars: 
horizontal-branch, stars: mass loss, stars: variables: other. Galaxy: formation, (Galaxy:) globular 
clusters: general, (Galaxy:) globular clusters: individual (M3, M13, NGC 1851), galaxies: dwarf 



1. Introduction 

Horizontal branch (HB) stars have long played a central role in the age debate. These 
low-mass stars, which burn helium in the core and hydrogen in a shell, are the immediate 
progeny of the luminous, vigorously mass-losing red giant branch (RGB) stars. Most 
importantly in the present context, their temperatures depend strongly on their total 
masses. More specifically, the lower the mass of an HB star, the bluer it becomes, by 
the time it reaches the zero-age HB (ZAHB). Therefore, the HB morphology in globular 
clusters (GC's) is naturally expected to become bluer with age. 

It has long been known that the first parameter controlling HB morphology is actually 
metallicity, with more metal-rich GC's present i ng re dder HB's than their more metal- 
poor counterparts. Still, ISandage fc WallersteinI (|l96C)[ ) first realized, based mainly on the 



ea rly observ ations of the GC's M3 ( NGC 5272) M13 (NGC 6205) , and M22 (NGC 6656) 
bv lSandagd (|l953l ). lArp fc Johnson! ()l955l) . and Urp fc Melbournd ([l95i), that GC's with 



a given metallicity might also present widely different HB types, due to the action of an 
unknown "second parameter." We quote from their study: 

"... the character of the horizontal branch is spoiled by the two clusters M13 and 
M22. (...) M13 appear[s] to be metal-rich, whereas the character of the horizon- 
tal branch simulates that of the very weak-lined group (M15, M92, NGC 5897). 
(...) M13 is younger than M2 or M5 (...) Consequently, in addition to chemical 
composition, the second parameter of age may be affecting the correlations." 

(Note that the sense o f the correlation between age and HB morphology suggested by 



Sandage fc WallersteinI is the opposite of what modern studies indicate to be necessary 
to account for the second-parameter phenomenon.) 

It soon became clear that age was not the only possible second parameter candidate. 
By the early 1970's, the list of candidates had increased sharply, and already inclu ded, in 



addit ion to age, the helium abundance and the abundances of the CNO elements ([Rood 



1973[ ). While the "age as the second pa rameter" scenario was to gain an important boost 



with the work bv lSearle fc ZinnI (|l978h . who noted that HB morphology tends to become 
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Figure 1. Comparison between differ ent mass los s recipes and the empirical d ata. The 
small crosses represent data from Judge & Stencel (1991,), sele cted according to [ Catelanl 
II2OOOI'). Fil led diamonds correspond to the HST-GHRS resuhs by iRobinson et alj 11998') and 
iMullan e Tal. (1998). Gray symbols with error bars represent data from the rec ent comp ilation 
bvlSchrSdcr & Cuntz (2007). In all panels, the solid lines show the fits derived bv'Cat elanI (i2(j60f ) 
from the Judge fc Stenceli data (using the different combinations of physical parameters indi- 
cated in the x-axis of each plot as the independent variable), except for panel e, where the fit is 
presented here for the first tim e . In p anel a, the dotted line represents the predicted mass loss 
rates according to the iReimerj (1197711 formula, w hereas the dashed line indicates the predicted 
mass loss rates according to lOrig lia et aT? (|2007|). In pane l e, the dasfi-dotted line indicates the 
mass loss rates predicted by the [Schroder fc Cunt j (|2005r i formula. 

redder with increasing Galactocentric distance - which was interpreted as an age effect, 
with more distant clusters being younger on average, and having possibly been accreted 
from "protogalactic fragments" of external origin over the Galaxy's lifetime - many other 
second parameter candidates have also surfaced over the years. This includes, among 
others, cluster concentration, total mass, and elli pticity; stellar rotation; magnetic fields; 
planetary systems; and mass loss on the RGB (seelCatelan''2008', for extensive references). 

While it seems clear now that age does play an important role, it has also become evi- 
dent that it is not the only parameter involved. Indeed, the presence of bimodal HB's in 
such GC's as NGC 2808 have l ong pointed to the need for other second parameters in ad- 
dition to age (e.g.. Rood et"al][l993 ). Recent, deep CMD studies have revealed that some 
of the most massive globulars, NGC 2808 included, present a surprisingly complex history 
of star formation, with the presenc e of extreme levels of helium enhanc ement among at 
least some of their stars fe.g.. iNorrij|2004HD'Antona et al.ll2005HPiotto et al..2005... 2007h . 
As noted by these authors, such He enhancement would provide a natural explanation 
for the presence of hot HB stars in these clusters. High helium abundances also appear 
to provide a natural explanation for some of the peculiarities observ ed in the CMD's and 
RR Lyrae properties in the GC's NGC 6388 and NGC 6441 fe.g.. ICatelan et al.ll2006l : 



Caloi fc D'Antonall2007l . and references therein). 

In the next few sections, we will address some empirical constraints that may be posed 
on some of these second parameter candidates. 
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Figure 2. Dependence of the integrated RGB mass loss upon age for fixed metallicity (panel a: 
[Fe/H] = -2.01; panel b: [Fe/H] = -1.31; panel c: [Fe/H] = -0.71) and upon metallicity for a 
fixed age (panel d: 9 Gyr; panel e: 12 Gyr; panel /: 15 Gyr), for the different mass loss recipes 
indicated. In panels a through c, the total mass loss has been normalized to a value of 0.20 Mq 
at 12 Gyr; in panels d through /, in turn, the integrated mass loss has been normalized to a 
value of 0.10 Mq at [Fe/H] = -1.54. 



2. Mass Loss in Red Giants 

In order to reliably predict the temperature of an HB star of a given composition and 
age, we must know how much mass it loses on the RGB. Unfortu nately, our know ledge 
of RGB mass loss remains rather limited. Most studies adopt the iReimersI (|l977[ ) mass 
loss formula to pre dict the i ntegrated mass loss along the RGB. However, recent evidence 
indicates that the iReimerg] formula is not a rehabl e description of RGB mass loss (e.g., 
Catelan|[2000l . l2008HSchroder fc Cu^l2005ll2007t ). In addition, there are several alter- 
native mass loss formulations which may better describe the available data. We illustrate 
this point by comparing, in Figure [U some empirical mass l oss rates with the predicted 
rates from several of these alternative mass loss formulae (see Catelan 2Q00l for extensive 
references). While the ,Reimers formula is clearly inconsistent with the data, the empirical 
data cannot conclusively distinguish among these alternative formulations 

The serious problem which this uncertainty in the RGB mass loss poses for our un- 
derstanding of HB morphology and its dependence on age and metallicity is apparent 



f Note that, while the lOriglia et al.] l|2007l ) mass loss formula is inconsistent with the plotted 
data, this should not be taken as evidence against its validity, since this formula has been sug- 
gested to apply exclusively to low-me tallicity stars. Still , some caveats regarding the lOriglia et al.l 
study have recently been raised (see lBover et al.ll2008D . 
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Figure 3. Comparis on between IStromgrenI l|l963f) pliotometry for M3, as derived by 
iGrundalil et al.l (|1998l . [1999), with the predictions of theoretical models for "Kms = 0.23, 
Z = 0.002. The observed data were shifted vertically so as to lead to a good match with 
the theoretical ZAHB at the red HB. 



from Figured! where the integrated RGB mass loss is plotted as a function of the age 
for fixed metallicity (panels a through c) and as a function of metallicity for fixed age 
(panels d through f). More specifically, we know that only a very mild AAf — [Fe/H] 
dependence can account for the observed relation between HB type and [Fe/H] with out 



resorting to a significant ag e-metallicity relation (se e, e.g., Fig, la inlLee et al. 1994 ). In 
this sense, we find that the Origlia et al.l ( 2007 ) and Schroder fc Cunt j l 2005 ) mass loss 
formulae lead to the weakest AM — [Fe/H] dependence. All other formu lae that we have 
tested lead to steeper dependencies between AM and [Fe/H] than the iReimer i (Il977l ) 
relation, thus implying steeper dependencies between age and metallicity as well. The 
precise dependence between AM and metallicity is also important in terms of explaining 
the u ltraviolet upturn phenomenon of elliptical galaxies and spiral bulges (see iCatelan 
I2OO7I for a recent review). 

Similarly, a stronger dependence between A A/ and age at fixed [Fe/H] makes it easier 
to account for a given second-parameter pair in terms solely of an age difference. Ac- 
cording to Figure [2l the equation that is most successful in this regard is a modified 



104 



M. Catelan 




Figure 4. Comparison between spectroscopically derived gravities for blue HB stars in M3 and 
M13 (from iBeiiil 12003 1 with tiie same theoretical models as before, but for two different Yms 
values: 23% (panel a) and 28% (panel 6). 



version of the ISchroder fc Cuntj (|2005l ) formula^ in which the adopted power law expo- 
ne nts are obtaine d by a least-squares fit to the Judge fc Stencel ( 199ll) data, selected as 
in Catelan ( 2000l ) (see Fig. [T]s). While not the steepest, the original ISchroder fc Cuntz 



formu la provides a stronger dependence between AM and age than does the iReimers 
(|l977f ) formula, which sho uld reduce the requir ed age difference between second parame- 
ter pairs. By contrast, the lOrigha et al.l (|2007l ) equation shows the weakest dependence, 
with a remarkably constant integrated AM value over a wide range in ages. 



3. Helium Enrichment in Globular Clusters 

As previously noted (§1), high levels of helium enrichment have been detected among 
some of the most massive Galactic GC's. Very recently, it has been suggested that such 
heliu m enhancements are in fa ct not the exception, but indeed the rule, among Galactic 
GC's (jP'Antona fc Caloill2008^ . Here we provide a first test of this scenario, in the case 
of the GC's M3 and NGC 1851. 

3.1. The Case of MS 

D'Antona fc Caloil (|2008l) and lCaloi fc D'Antonal (|2008l ) claim that the blue HB compo- 
nent in MS owes its origin to a moderate level of He enhancement in the cluster, between 
2% and 6%. Is this supported by the available data? 

To answer this q u estion , we compare, in Figure 5 cano nical theoretical predictions 
from ICatelan et"aLl (|l998[) and ISweigart fc Catelann|l998[ ) for a helium abundance of 
Ym.s = 0.23 and a metallicity Z = 0.002 wit h high-precision photometry, in the'Strom grenI 



19631 ) system, fro miGrundahl et al. ( 19981 [1999 ). The empirical data were corrected for 



eddening following Harris! ( 1996 ). In these plots, the drawn lines represent different fidu 



cial loci, as follows: ZAHB, middle-age HB (MAHB, or average locus occupied by the 
HB stars), 90%-age HB (90AHB, or locus below which one should expect to find about 
90% of all HB stars), and terminal-age HB (TAHB, or He exhaustion locus). Except 
for a discrepancy between the predicted and observed numbers of highly evolved stars 
(both on the blue and red HB sides of the distribution) that was previously noted by 
Valcarce fc CatelanI (|2008l ). one finds remarkable agreement between the model predic- 
tions for a constant Y a nd the observations. Such an agreement is also confirmed by the 
spectroscopic data from Behij ((2003) , as can be seen from Figure SI where we limit the 
comparison to temperatures lower than 11,500 K due to the well-known complications 
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Figure 5. The age difference that is required to account for the difference in HB morphology 
between the blue and red components of NGC 1851's HB {lines, based on the different mass 
loss formulae discussed in 321) is compared with the age difference that is estimated from the 
observed split on the SGB {gray band). 



brought about by the "Grundahl jump" phenomenon (jGrundahl et al. Interest 



ingiy, this plot also appears to support a similar helium abundance between M3's blue 
HB stars and the redder blue HB stars in M13. 

3.2. The Case of NGC 1851 



Milone et al.l (j2008f ) have recently discovered that the subgiant branch (SGB) of NGC 1851 
is actually split into two separate components, which may be linked to the cluster's well- 
known bimodal HB morphology. The most straightforward explanation for this split 
would be a difference in age by 1.0 ± 0.4 Gyr. However, as shown in Figure [5l this is 
inconsistent with the age difference that would be required to fully account for the sepa- 
ration between the blue and red HB components of the cluster, irrespective of the mass 
loss formula (SQI) used. A difference i n metallicity between th e two components is also 



ruled out by recent spectroscopic data ( Yong k, Grundahll2008l) . Here we apply the same 



CMD test as in the previous section to constrain the possibility of a difference in Y be- 
tween the two components as being responsible for the well-known bimodal nature of the 
cluster's HB. 

The result is shown in Figure [6l While the quality of the data is not as high as in the 
case of M3, one is still able to derive some general conclusions. First, the same theoretical 
ZAHB does appear to provide a reasonable description of the lower boundary of the data, 
both for the red and blue HB components - which suggests that at least some of the stars 
on the blue HB have the same Y as do the red HB stars. Second, there is a predominance 
of ovcrluminous stars on the blue HB, at colors around {b—y)o « 0.05 — 0.15. While these 
might in principle be interpreted in terms of a moderate level of helium enrichment, 
perhaps of the order of 3% — 4% on average (see Fig. [T]), the more straightforward 
explanation is that these stars actually represent the well-evolved progeny of the blue 
ZAHB stars that are found at higher temperatures. If so, this would again suggest that 
the late stages of HB evolution are somehow significantly undere stimated by present- 



day HB tracks, similar t o what was previously found elsewhere (jCatelan et al.l 12001 
IValcarce fc Catelan 2008 . and references therein). 



Clearly, more work is needed before we ar e able to conclusive l y establish the nat ure of 
NGC 1851's bimodal HB and SGB (see also lCassisi et al.ll2008t Isalaris et al.ll2008l ). 
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4. The OosterhofF Dichotomy and the Formation of the Milky Way 

Irrespective of our ability to properly model HB stars, we can use RR Lyrae stars to 
derive entirely empir ical constraints on the process of formation of the Milky Way. In the 



Searle fc ZinnI (|l978f ) scenario, much like in modern ACDM cosmology, one expects galax- 



ies such as our own to have formed by the accretion of "protogalactic fragments" that 
may have resembled the early counterparts of the Milky Way's present-day dwarf satel- 
lite galaxies. Useful constraints on recent accretion events may be posed by the presence 
of younger stellar populations in several of these galaxies ( Unavane et al.lll996l ). Still, in 
order to probe what happened very early on, we must look a t the a ncient components - 
and RR Lyrae stars are especially useful in that regard fe.g.. ICatelan 2008 ). 

Arc the ancient populations in the Milky Way's dwarf satellites, as traced by their 
RR Lyrae pulsators, consistent with the Galactic spheroid having been built therefrom? 
The answer is provided in Figure [51 where we compare the average properties of the 
fundamental-mode (ab-type) RR Lyrae stars in Galactic {left panel) vs. nearby extra- 
galactic {right panel) GC's and field populations. While the Galactic distribution clearly 
presents the so-called Oosterhojf dichotomy, with a tendency for syste ms to clump around 
the "Oosterhoff I" (Ool) and "Oosterhoff 11" (OoII) regions (see also lMiceh et al.ll2008i 
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for the case of halo field stars), the opposite happens in the case of nearby extragalactic 
systems, which tend to be preferentially Oosterhoff-intermediate. This strongly suggests 
that the oldest components of the Galaxy cannot have been formed by accretion of even 
the early counterparts of its present-day dwarf galaxy satellites. 

As indicated in Figure [8l at l east one of the newly discovered S PSS dwarf galaxie s 
fe.g. jBelokurov et al.ll20 06n2007|) , CVn I, is Oosterhoff interinediat e (|Kuehn et alfeOOSi ). 
whereas the Bootes dwarf is OoII (iDairOra et al.ll2006l : [Siegell200d) . Unfortunately, some 
of the low-mass SDSS galaxies seem to harbor a mere one or two RR Lyrae stars, which 
makes it more difficult to assign a conclusive Oosterhoff status to them. Indeed, due to 
statistical fluctuations, and since the HB lifetime is of order 100 Myr, it is not entirely 
clear whether the same Oosterhoff types would necessarily be inferred for these galaxies 
if they were observed, say, a few hundred Myr in the future (which is very little, in terms 
of Galactic history), when these HB stars will have long left the HB phase, to be replaced 
by an entirely new generation of HB stars. Still, the present-day properties for several 
of the RR Lyrae stars that are fou nd in these very lo w-mass galaxies do appear to be 
consistent with an OoII status fe.g.. iGreco et al. I l2008h . 
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Figure 8. Distribution of Galactic GCs (left) and stellar populations associated with neigh- 
boring dwarf galaxies (right) in the average ab-type RR Lyrae period (Pab) vs. [Fe/H] plane. 
Ga lactic GC's are classified into "bu lge/d isk," "young ha lo," and "old halo" subsystems follow- 
ing [MickixASndinBiEES (HqqI). See lCatelanI (|2008l ) for further details. 



5. Conclusions 

• HB stars play a central role in the age debate. Still, before we are able to predict 
how (ZA)HB temperature changes with age, we must properly describe RGB mass loss. 

• In studies of HB morphology, it is not sufficient anymore to analyze solely the "hori- 
zontal" HB morphology, meaning the temperature and/or color distribution of HB stars: 
"vertical" HB morphology, or the distribution of HB stars in luminosity at fixed tem- 
perature (or color), provides us with unique information to help us constrain theoretical 
scenarios for the origin and evolution of these stars. Indeed, the available data appears to 
strongly constrain, if not conclusively rule out, the possibility of significant He enhance- 
ments among M3's blue HB stars, while at the same time suggesting that canonical HB 
models underestimate the duration of the late stages of HB evolution. 

• Irrespective of our ability to model them, HB stars - and RR Lyrae in particular - 
represent invaluable tools to probe into the Milky Way's early formation history. 

Acknovirledgments. I would like to warmly thank Gisella Clementini, Frank Grundahl, 
Bob Rood, Horace Smith, Allen Sweigart, and Aldo Valcarce for useful discussions, com- 
ments and suggestions. Support for this work is provided by Proyecto Basal PFB-06/2007, 
by FONDAP Ccntro de Astroffsica 15010003, by Proyecto FONDECYT #1071002, and 
by a John Simon Guggenheim Memorial Foundation Fellowship. 



References 

Arp, H. C., & Johnson, H. L. 1955, ApJ, 122, 171 
Arp, H. C., & Melbourne, W. G. 1959, AJ, 64, 28 
Behr, B. B. 2003, ApJS, 149, 67 
Belokurov, V., et al. 2006, ApJ, 647, Llll 
Belokurov, V., et al. 2007, ApJ, 654, 897 

Boyer, M. L., McDonald, I., van Loon, J. Th., Woodward, C. E., Gehrz, R. D., Evans, A., & 
Dupree, A. K. 2008, AJ, 135, 1395 



HB Stars and The Age Debate 



109 



Caloi, v., & D'Antona, F. 2007, A&A, 463, 949 
Caloi, v., & D'Antona, F. 2008, ApJ, 673, 847 

Cassisi, S., Salaris, M., Pietrinferni, A., Piotto, G., Milone, A. P., Bedin, L. R., & Anderson, J. 

2008, ApJ, 672, L115 
Catelan, M. 2000, ApJ, 531, 826 

Catelan, M. 2007, in: New Quests in Stellar Astrophysics II: The Ultraviolet Properties of 

Evolved Stellar Populations, in press (astro-ph/0708.2445) 
Catelan, M. 2008, to appear in Ap05S'(astro-ph/O5O7464) 

Catelan, M., Bellazzini, M., Landsman, W. B., Ferraro, F. R., Fusi Pecci, F., & Galleti, S. 2001, 
AJ, 122, 3171 

Catelan, M., Borissova, J., Sweigart, A. V., & Spassova, N. 1998, ApJ, 494, 265 

Catelan, M., Stetson, P. B., Pritzl, B. J., Smith, H. A., Kinemuchi, K., Layden, A. C, Sweigart, 

A. v., & Rich, R. M. 2006, ApJ, 651, L133 
DaU'Ora, M., et al. 2006, ApJ, 653, L109 

D'Antona, F., Bellazzini, M., Caloi, V., Fusi Pecci, F., Galleti, S., & Rood, R. T. 2005, ApJ, 
631, 868 

D'Antona, F., & Caloi, V. 2008, MNRAS, 390, 693 

Greco, C, DaU'Ora, M., Clementini, G., et al. 2008, ApJ, 675, L73 

Grundahl, F., Catelan, M., Landsman, W. B., Stetson, P. B., & Andersen, M. I. 1999, ApJ, 524, 
242 

Grundahl, F., VandenBerg, D. A., & Andersen, M. I. 1998, ApJ, 500, L179 

Harris, W. E. 1996, AJ, 112, 1487 

Judge, P. G., & Stencel, R. E. 1991, ApJ, 371, 357 

Kuehn, C, Kinemuchi, K., Ripepi, V., et al. 2008, ApJ, 674, L81 

Lee, Y.-W., Demarque, P., & Zinn, R. 1994, ApJ, 423, 248 

Mackey, A. D., & van den Bergh, S. 2005, MNRAS, 360, 631 

Miceh, A., Rest, A., Stubbs, C. W., et al. 2008, ApJ, 678, 865 

Milone, A. P., Bedin, L. R., Piotto, G., & Anderson, J. 2008, A&A, in press (astro-ph/0810.2558) 
MuUan, D. J., Carpenter, K. G., & Robinson, R. D. 1998, ApJ, 495, 927 

Munoz, R. R., Carhn, J. L., Frinchaboy, P. M., Nidever, D. L., Majewski, S. R., & Patterson, 

R. J. 2006, ApJ, 650, L51 
Norris, J. E. 2004, ApJ, 612, L25 

Origlia, L, Rood, R. T., Fabbri, S., Ferraro, F. R., Fusi Pecci, F., & Rich, R. M. 2007, ApJ, 667, 
L85 

Piotto, G., Bedin, L. R., Anderson, J., et al. 2007, ApJ, 661, L53 
Piotto, G., Villanova, S., Bedin, L. R., et al. 2005, ApJ, 621, 777 
Reimers, D. 1977, A&A, 57, 395 

Robinson, R. D., Carpenter, K. G., & Brown, A. 1998, ApJ, 503, 396 
Rood, R. T. 1973, ApJ, 184, 815 

Rood, R. T., Crocker, D. A., Fusi Pecci, F., & Buonanno, R. 1993, in: G. H. Smith & J. P. 

Brodie (eds.), Proc. ASP Conf. Ser., Vol. 48, The Globular Cluster- Galaxy Connection, 

(San Francisco: ASP), p. 218 
Salaris, M., Cassisi, S., & Pietrinferni, A. 2008, ApJ, 678, L25 
Sandage, A. R. 1953, AJ, 58, 61 
Sandage, A., & Wallerstein, G. 1960, ApJ, 131, 598 
Schroder, K.-P., & Cuntz, M. 2005, ApJ, 630, L73 
Schroder, K.-P., & Cuntz, M. 2007, A&A, 465, 593 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 
Siegel, M. H. 2006, ApJ, 649, L83 
Stromgren, B. 1963, QJRAS, 4, 8 
Sweigart, A. V., & Catelan, M. 1998, ApJ, 501, L63 
Unavane, M., Wyse, R. F. G., & Gilmore, G. 1996, MNRAS, 278, 727 
Valcarce, A. A. R., & Catelan, M. 2008, A&A, 487, 185 
Yong, D., & Grundahl, F. 2008, ApJ, 672, L29 



